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ABSTRACT 



We report results from a multi-wavelength analysis of the Galactic SNR RCW 86, 
the proposed remnant of the supernova of 185 A.D. We report new infrared observa- 
tions from Spitzer and WISE, where the entire shell is detected at 24 and 22 /mi. We 
fit the infrared flux ratios with models of colli sionally heated ambient dust, finding 
post-shock gas densities in the non-radiative shocks of 2.4 and 2.0 cm" 3 in the SW 
and NW portions of the remnant, respectively. The Balmer-dominated shocks around 
the periphery of the shell, large amount of iron in the X-ray emitting ejecta, and lack 
of a compact remnant support a Type la origin for this remnant. From hydrodynamic 
simulations, the observed characteristics of RCW 86 are successfully reproduced by 
an off-center explosion in a low-density cavity carved by the progenitor system. This 
would make RCW 86 the first known case of a Type la supernova in a wind-blown 
bubble. The fast shocks (> 3000 km s" 1 ) observed in the NE are propagating in the 
low-density bubble, where the shock is just beginning to encounter the shell, while 
the slower shocks elsewhere have already encountered the bubble wall. The diffuse 
nature of the synchrotron emission in the SW and NW is due to electrons that were 
accelerated early in the lifetime of the remnant, when the shock was still in the bubble. 
Electrons in a bubble could produce gamma-rays by inverse-Compton scattering. The 
wind-blown bubble scenario requires a single-degenerate progenitor, which should 
leave behind a companion star. 



1. Introduction 



It is widely believed that Type la supernovae (SNe) originate from white dwarf (WD) 
stars that have been pushed close to the Chandrasekhar limit, either by accreting matter from a 
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companio n star (sing 



scenario) (|lsern et al 



e-deg enerate scenario), or by merging with another WD (double-degenerate 



2008). The detailed physics of the single-degenerate scenario, particularly 



the role, if any, that the binary s ystem plays in shaping the surrounding medium is poorly 



understood. 



Hachisu et al. 



(1996) suggested that, under certain conditions, accreting WDs in 
a binary system with a main sequence or red giant companion might blow a substantial wind, 
carving out low-dens ity cavities in the surrounding ISM in the ~ 10 6 yr prior to explosion. 



Badenes et al 



(|2007|) considered young (< 1000 yr) supernova remnants (SNRs) known to be of 



Type la origin, finding that none were consistent with evolution into such a modified medium. 

The Galactic SNR RCW 86 (G315.4-2.3; MSH 14-63) is one of a small number of 
remnants with a historical c onnection, as it has been proposed to be the remains of SN 185 A .D. 



(Cla rk & Stephenson 



1975|). The distance has been reporte d to be near 1 kpc (IKaastra et al. 



992), but more re cent measurements have put it from 2.3 (|Sollerman et al 



(IRosado et al 



2003b to 2.8 kpc 



1996). Throughout this work, we adopt a distance of 2.5 kpc, para meterized, wher e 



appro priate, as di$- The claim that this event was a SN at all has been disputed (|Chin & Huang 



1994), but 



Stephenson & Green 



(|2002|) review the entirety of the available information on the 



g uest star" o f 185 A .D., concluding that it was indeed a SN, a view supported by the work 



of 



Zhou et al. 



(20061), a nd th at RCW 86 represents the most likely candidate for its remnant. 



According to iSchaeferi (| 19931) . the maximum unrefracted altitude of RCW 86 as seen from the 
city of Luoyang, China, was 2.3°, leading to an airmass of 18. Under good seeing conditions 
(0.2 magnitudes extinction per airmass), the atmospheric exti nction would be A v = 3.6. O ptical 



measurements have placed the ISM extinction at ~ Ay = 1.7 (|Leibowitz & Danzigei 



1983b, so 



a supernova at 2.5 kpc (see below) with, for instance, M v = -17 would appear as a new star of 
m = 0.3, easily visible to ancient observers. Significant sources of error are possible in virtually 
all of these parameters. If the Chinese historical observations came from a region farther south 
in China, the airmass in the direction of RCW 86 would be significantly reduced . For instance, 



an altitude of 5° above the horizon implies 10 airmasses (IKasten & Young 



19891) . reducing the 
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atmospheric extinction to 2 mag. A brighter Type la SN, which we argue in favor of later, could 
add another 2 mag to the observed brightness. On the other hand, we report X-ray column 
densities of ~ 6 x 10 21 cm" 2 , which could raise the interstellar extinction from 2 to 4 mag. In 
summary, it appears quite possible that there could have been a supernova visible from China at 
the position of RCW 86, at a distance of 2.5 kpc. 

The recent measurement of high proper motion of the shock, as defined by the on set of X-ray 



emission in the NE limb implies very fast shocks (~ 6000 km s ! ) (|Helder et al 



2009 ), supporting 



a young age for the remnant. The broad Fe Ka line width me asured by Suzaku (lUeno et al 



and the presence of non -thermal synchrotron X-ray emission (Ba mba et al. 



2001a 



Vink et al. 



2000; 



2007) 



Borkowski et al 



2006|) als o suggest a young SNR. The remnant shows a complete she 



optical, and X-ray energies (iKesteven & Caswel 



19871 : 



Smith 



1997; 



Pisarski et al. 



1 at radio, 



1984), and has 



20091). Shock speeds 



recently been detected in TeV gamma-rays by H.E.S.S. (|Aharonian et al.1 
vary greatly from one side of the remnant to the other, from the high shock speeds in the east 
to the slower (550 - 650 km s" 1 ) shocks in the NW and SW corners, as de ri ved from opt i cal 



spectroscopy (|Long & Blair 



1990 



Ghavamian 



1999 



Ghavamian et al. 



2001). 



Helder et al. 



(2009) 



also report that the eastern limb could be a site of significant cosmic-ray acceleration, based on 
the discrepancy between proper motions of the X-ray emitting filaments and the spectroscopically 
determined proton temperature. 

Numerous authors have pointed out the difficulties with reconciling the young age of RCW 
86 with its large size. At a distance of 2.5 kpc, the average speed of the shock for this large 
remnant (40' in diameter) is ~ 7800 km s" 1 , but as previously mentioned, shock speeds in most 
of the remnant are an order of magnitude less than this. It has been suggested (e.g., Vink et 
al. 1997) that a cavity explosion offers the most natural explanation, with the shock racing 
through a low-density bubble, recently having encountered a dense shell. Cavity explosions are 
typically associated with core-collapse (CC) supernovae, but the nature of the progenitor system 
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for RCW 86 is uncertain. iBadenes et al.1 (I2007|) suggested that RCW 86 may be the result of a 
Type la explosion into a cavity formed by the progenitor system, a scenario that could offer new 
constraints on accretion wind outflows from single-degenerate progenitor systems of Type la SNe. 

This paper is organized as follows: In Sections 2 and 3, we report the results of new infrared 
(IR) observations with Spitzer and WISE. We discuss the nature of the progenitor system in 
Section HI and the evidence for a cavity explosion in Section |5] The nonthermal synchrotron 
emission is discussed in Section 6, and the gamma-ray emission in Section 7. 



2. Infrared Observations 



IR emission from SNRs can be thermal continuum from warm dust grains heated by 
collisions with energetic ions and electrons in the X-ray emitting gas, line emission from low to 
medium-ionization states of atoms cooling in the post-shock region, or a combination of the two. 
Nonthermal emission via synch rotron radiation c an contribute as well, but is generally negligible 



compared to thermal processes (IRho et al. 



2008). Generally speaking, emission from younger 



SNRs expanding into a low-densi ty medium (characte rized by non-radiative shocks in the optical) 



will be dominated by continuum (IWilliams et al. 



201 1|) . while older shocks or th ose encounte ring 



molecular clouds will become radiative and show strong mid-IR line emission (IHewitt et al. 



20091). IR studies of SNRs offer the possibility to learn about both the conditions of the remnant 



and the physics of the gas-dust interaction in shocks. 



Previous IR observations of RCW 86 have led to conflicting results. iDwek et al 



(1987) used 



IRAS observations to determine that the IR/X-ray flux ratio (IRX ratio) from the bright SW shell 
(the only part of the remnant firmly detected by IRAS) is higher than theoretical predictions. They 
attributed this to either a dust/gas mass ratio higher than average for t he Galaxy or a poor coupli ng 



between the IR emitting dust and the X-ray emitting gas. However, 



Greidanus & Stroml ( 1990) 
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used the same IRAS data and IR emission model to arrive at the opposite conclusion, namely that 
the IRX ratio for the remnant is lower than predictions by a factor of 10, implying a low dust/gas 
mass ratio. 

Studies done with IRAS suffer from both poor angular resolution and low sensitivity, and 
resolving discrepant results such as these requires study by advanced IR telescopes. Using Spitzer 
and WISE, we detect the entire shell at 24 and 22 //m and the SW and NW portions at 12 and 70 
fim. Faint extended emission in the SW is seen in short-wavelength bands as well. Detection of 
the eastern shell of the remnant is significantly complicated at 70 /im by confusion with Galactic 
IR emission. 

As part of Spitzer program 50698, we obtained a full map of the remnant at 24 and 70 /im 
using the Multiband Imaging Photometer for Spitzer (MIPS) on 2009 March 3 1 , just weeks before 
the end of the cold Spitzer mission. The size of the remnant (~ 40' in diameter) required two 
scan maps (AOR ID 26750208 & 26749952) of 1° x 0.5°, which we mosaicked together using 
the MOPEX software provided by the Spitzer Science Center. We used a medium scan rate for 
the maps, offsetting each scan leg by 148" to ensure full coverage at 70 /im and redundancy by a 
factor of 2 at 24 pm. We used all four channels (from 3.6 to 8 pm) of the Infrared Array Camera 
(IRAC) to map a region of ~ 7' x 7' on 2009 March 23 (AOR ID 26748928), centered on the 
brightest filaments in the SW region of the shell seen in X-rays with Chandra (Rho et al. 2002, 
hereafter R02). Both the IRAC and MIPS data were processed using version 18.7 of the Spitzer 
pipeline. The pipeline processed 70 fim data contains some artifacts, most notably striping in the 
direction of the scan leg. 

In mid- April of 201 1 , the Wide-Field Infrared Survey Explorer (WISE) science team released 
the first half of the data taken by WISE, a mid-IR telescope which surveyed the entire sky at 3.3, 
4.6, 12, and 22 pm. RCW 86 was in the region of the sky covered by this initial data release. To 
provide a more complete picture of the IR properties of the SNR, we have downloaded the WISE 
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images of this remnant at all four wavelengths. 



3. Analysis and Modeling 
3.1. Morphology 

We show the MIPS images at 24 and 70 /im in Figured] along with X-ray and optical images 
as described in the caption. At 24 /mi, where the angular resolution of the telescope is ~ 7", the 
correlation between the IR, X-ray, and optical images is significant, implying a tight coupling 
between the IR emitting dust and the X-ray emitting gas. A strong correlation exists for the 70 /im 
data (resolution ~ 20") as well in the SW and NW portions of the shell; however, the eastern part 
of the shell is dominated by emission that is likely from a foreground or background region. The 
IR colors in the E region allow us to trace the forward shock and make a determination of what 
is and is not associated with the remnant. At 24 /im, the much lower background level allows 
distinction between this region and the thin, non-radiative filaments of the blast wave. 

We show a three-color IRAC image in Figure [2] where a few faint filaments are observed 
that correspond to fe atures seen at 24 /mi. These filaments also appear as radiative filaments 



in the [S II] image (|Smith 



1997 ). making it likely that the IRAC emission comes from lines 
(most likely [Ar II] at 7 /im and [Fe II] at 5.4 /im). It is possible that emission from polycyclic 
aromatic hydrocarbons (PAHs) heated immediately behind the shock front could contribute as 
well, although the lifetime of these molecules is quite short in the post shock environment. 

A mosaic of the WISE 4.6 /mi, 12 /mi, and Spitzer 24 /im images is shown in Figure[3] WISE 
data show the entire periphery of the shell at 22 /mi (though with a lower angular resolution than 
with Spitzer), as well as significant portions of the SW shell at 12 /im, where the WISE resolution 
is comparable with Spitzer's at 24 /im. Faint emission from the NW filament is detectable at 12 
/mi as well. At 4.6 and 3.3 /im, only a small filament is visible in the extreme SW, a filament 
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which was not covered by the spatially limited IRAC map. The colors of this image, which reflect 
the temperature of the dust, separate the emission that is associated with the SNR from that of the 
background. This is particularly important in the east, where the thin red filaments trace the blast 
wave. 

3.2. IR Fluxes and Ratios 

For our analysis, we selected two regions in the remnant with detections at both 24 and 70 
/mi: a small region in the bright SW "corner" that corresponds to purely non-radiative emission, 
as identified in optical images, and a filament in the NW stretching ~ V. The regions are indicated 
in Figure Q] and Figure [2l and are significantly larger than the angular resolution of Spitzer at 
70 fim. To measure the flux, we chose several local background regions outside of the remnant 
near the region of interest. The background is somewhat variable, particularly in the NW, so we 
average the background subtracted fluxes and report the standard deviation of these measurements 
as the error. Because we are interested in the emission from dust only, we must account for both 
stars and line emission from cooling gas. Neither region contains bright stellar sources, and only 
the 24 fim image shows any stars at all in the region of interest. We excise regions at 24 /im 
where stars are seen, using IRAC (where available) and optical images to determine the stellar 
nature of unresolved sources. Despite no visible stars in the regions at 70 fim, we excise the same 
regions when measuring the 70 fim flux. We find that stellar contamination to the 24 /mi flux is 
negligible in the small SW region and accounts for 15% in the larger NW region. The integrated 
flux from the excised regions at 70 fim was less than 1% of the total. Fluxes reported in this work 
are the "star-subtracted" fluxes. We measure fluxes in the SW of 0.37 and 1.41 Jy at 24 and 70 
fim, respectively, and 0.69 and 3.10 Jy in the NW. Errors on these flux measurements are ~ 6% 
in the SW and ~ 30% in the NW, where the larger errors in the NW result from a fainter, more 
diffuse source region with a more variable background. The 70/24 fim flux ratios are 3.8±0.3 and 
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4.5±2 in the SW and NW, respectively. 

To account for the contribution of line emission to the integrated fluxes, we use IRS spectra 
taken as part of our Spitzer observing program. We defer analysis of the IRS data to a future 
publication, but note that a spatially integrated low-resolution spectrum from the SW shell 
shows the presence of the [Fe II]/[0 IV] line blend at ~ 26 /mi (the lines are unresolved by the 
spectrograph, which has spectral resolution A/AA ~ 100). This line lies on top of a strong dust 
continuum, and when we integrate the IRS spectrum over the MIPS 24 bandpass and response 
curve, the total line flux accounts for 1.2% of the observed 24 /mi flux. The extraction region 
for this spectrum is shown in the bottom left panel of Figured! and lies on top of both radiative 
and non-radiative filaments. Although we do not have IRS spectra of the NW filament, radiative 
shocks there are fainter, and we expect that line contamination to the 24 /mi fluxes should be even 
less important than in the SW. We do not have spectral data extending to the wavelengths of the 
70 /mi detector, but the bandpass of that filter is extremely broad (FWHM ~ 20 /mi), and any 
lines present, such as [O I] at 63 fim or [O III] at 88 fim (where filter throughput is only 14%) 
would have to be extraordinarily strong to contribute significant flux. We thus ignore any line 
contributions to the fluxes reported above. 



3.3. Modeling 



To model IR emis s ion, we follow a pr ocedure identical to that described in previous works 



( B orko wskietal. 



2006 



Williams et al. 



2006J). The heating rate of a grain immersed in a hot 



plasma depends on both the post-shock gas density (collision rate) and temperature (energy per 
collision). In general, the temperature of a plasma can be determined by fits to X-ra y spectra. 



Dust heating models depend on this temperature, but this dependency is not strong (jDwek 



1987) 



The density of the gas is a more sensitive parameter for dust heating models and controls the 
temperature of the grains (i.e. the 70/24 /mi flux ratio). We regard the gas density as a free model 
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parameter and adjust it to match the observed IR fluxes. 



We account for sputtering of grains by ions, using a plane-shock model that superimposes 

ft 

regions of increasing sputtering timescale, defined as r sp = J n p dt, where n p is the post-shock 
proton density, while keeping temperature behind the shock constant. (Note that the sputtering 
timescale is directly related to the ionization timescale used in X-rays, defined as t, = J^n e dt.) 
Inputs to the model are a pre-shock grain size distribution, grain type a nd abu ndance (for these 



quantities, we use models for the Galaxy from 



Weingartner & Draind (|200l|) . which contain 



a mixture of amorphous silicate and graphite grains), and ion and electron temperature and 
density. The code calculates heating and sputtering for grains from 1 nm to 1 /mi, producing a 
unique temperature and spectrum for each grain size. The spectra are summed according to the 
post-sputtering size d istribution to produce a single model spec trum. We use sput tering rates from 



Nozawa et al. 



Jurac etal 



d 19981) . Although we 



(|2006|) with enhancements for small grains from 
do not include it in our calculations, w e note that the shock speeds in this remnant are only slightly 

(2010) as 500 km s" 1 , where nonthermal, or 



above the upper threshold, reported by 



Sankrit et al 



ballistic sputtering effects begin to become important. 

To determine kT and r, for the X-ray emitting gas, we analyze X-ray spectra from archival 
Chandra and XMM-Newton observations. The Chandra data were obtained in 2001 (Obs. ID 
1993), where the SW shell was observed for 92 ks. We reprocessed and extracted the data using 
version 4.3 of the Chandra Interactive Analysis of Observations (CIAO). We obtained archival 
X-ray data in the NW rim and SW corner from the XMM-Newton EPIC-MOS detectors, observed 
in 2007 (Obs. ID 0504810401 & 0504810301, J. Vink, P.I.). The XMM-Newton data were 
processed using version 9.0 of the Science Analysis Software (SAS) for XMM. We filtered the 
data from MOS 1 and MOS 2 for flares and bad pixels, group ing the data to a minimum of 25 



counts per bin. We fit the data using version 12.6.0 of XSPEC (jArnaud 



1996). 
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3.4. SW 



R02 used the high spatial resolution of Chandra to uncover a complex system of shocks in 
the SW, with hard regions from the shocked ejecta mixed with softer regions from shocked ISM 
and regions dominated by non-thermal emission. We employ a similar modeling procedure to that 
used in R02 to fit the X-ray spectra from our region in the SW. We jointly fit Chandra ACIS and 
XMM-Newton EPIC-MOS 2 data (this particular region falls on the dead MOS 1 chip), letting the 
abundances of C, N (tied to C), O, Mg, Si, and Fe float freely. We find a reasonable fit with kT 
= 0.39 keV and r, = 6.0 x 10 10 cm" 3 s. The reduced x 2 value in this simple model is rather high, 
~ 3. R02 found similarly poor values of x 2 m modeling emission from this region, attributing it 
to poor Fe L-shell data in low temperature plasmas (see discussion in Section 3.2 of that paper). 
Excising the regions of the spectra that were poorly fit improved the value of x 2 , but did not 
affect the fitted parameters of the thermal model, which are fully listed in Tabled] Freezing these 
parameters (T e and r,) in our dust heating model, we find that the IR luminosity and flux ratio are 
best fit with a post-shock density of n H = 2.4 (2.1, 2.75) cm" 3 and a radiating dust mass of 3.2 d\ 5 
(2.7, 3.9) x 10" 4 M Q . Lower and upper limits, listed in parentheses, are derived from fitting dust 
emission models to the extrema of the 5 70 /5 2 4 ratio. We assume that T p = T e = 0.39 keV, which i s 



appropriate for the shock speeds and ages seen in this section of the remnant (|Ghavamian et al. 



2007 



van Adelsberg et al. 



entirely by electrons (|Dwek 



2008). In any case, grain heating at these temperatures is done almost 



1987) 



3.5. NW 



We show the 0.4-7 keV EPIC-MOS spectra from the NW in Figure H The shock speed 

I 1 1 I 

of this filament is 580-660 km s" 1 (|Ghavamianl 1 1 999|) . Visible in the spectra is a continuum 
extending to high energies, featureless except for an Fe Ka line a t ~ 6.4 keV. We account fo r these 



components with an srcut model for the synchrotron emission (|Reynolds & Keohane 



1999) and 
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an Fe-only vpshock model with kT = 5 keV and r, = 10 9 cm" 3 s. R02 also found Fe Ka emission 
in the SW, where shock speeds and plasma temperatures are also relatively low, attributing it to 
reverse-shocked ejecta. It is interesting to note that the nonthermal continuum in the NW does not 
show the same morphology as the low-energy thermal emission. Thermal X-rays in that region are 
well contained to the thin filament seen in Figured] but the synchrotron emission appears much 
more diffuse. We show this in Figure [51 Similar discrepancies were observed in the SW, leading 
R02 to conclude that the synchrotron emission is associated with reverse shocks driven into the 
ejecta, and not the forward shock. 

We follow the same modeling procedure as for the SW, and fit the spectra with a vpshock 
model (with high-energy components added as described above). In this region, we jointly fit 
MOS 1 and 2 data, as this part of the remnant has not been observed with Chandra. The best fit 
model (reduced x 2 = 2.20) has kT = 0.28 (0.27, 0.30) keV and r ; - = 9.9 (8.7, 1 1.1) x 10 10 cm" 3 s, 
with an emission measure (EM) for the thermal component of 2.51 d\ 5 (1.90, 2.91) x 10 57 cm" 3 
(2.15 dl 5 (1.63, 2.49) cm" 3 M ). We let the abundances of C, N (tied to C), O, Mg, Si, and Fe 
float freely in the model, and find that all are in the neighborhood of solar (using 



Wilms et al. 



(2000) abundances), with the exception of Si (1.7 solar) and Mg (0.6 solar). Ne was kept fixed to 
solar as an anchor for the other abundances, as Ne is not expected in the ejecta of a Type la SN, 
and is not depleted onto grains in the ISM. The overabundance of Si is likely a result of using 
a single-temperature shock model to fit the data. As we show in Section [51 the shock structure 
is likely more complicated than this, with varying shock speeds, and a single-temperature shock 
model with kT = 0.28 keV could underpredict the Si Ka line at 1.8 keV. The model compensates 
for this by making Si overabundant. Mg underabundance might be caused by depeletion onto dust 
grains. 

We find a post-shock density from matching IR data of n H = 2.0 (0.75, 5.0) cm" 3 and a 
radiating dust mass of 1.2 d\ 5 (0.33, 3.6) x 10" 3 M . The large uncertainties for this region are a 
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result of the errors in the flux measurement (see Section 13. 21) . Th e plasma temperature derived 



from this model is significantly lower than that of 



Vink et al. 



(|2006h . who fit XMM spectra from 



a similar region in the NW with a two-component NEI model , finding temperatu res of ~ 1 and 



Helder et al. 



(1201 II) . although they 



4 keV. However, our temperature is close to that reported by 
report lower ionization timescales. The difference is likely due to the fact that our "NW" spectrum 
is extracted from a region that is farther east than theirs, and much brighter in both X-rays and IR, 
implying a higher density. Differences also result from different spectral models used. 

It is worth noting that the ionization timescales obtained from fitting X-ray spectra, combined 
with the post-shock densities from fitting IR data, allow a determination of the age of the shock. 
We obtain ages of ~ 700 and ~ 1300 yr for the SW and NW, respectively. The se ages, derived 



from plane-shock model fits, are underestimates for the true age of the remnant ( Borkowski et al. 



2001b)), but do strongly favor a young SNR, strengthening the case for RCW 86 as the remnant of 
SN 185 A.D. Where dust is present and X-ray spectra of the swept-up ambient medium can be 
examined, this method provides a useful check on the ages of young SNRs. 

In modeling dust emission from the post-shock environment, we focus only on the Spitzer 
data, as the wavelengths covered with Spitzer are more appropriate for the bulk of the dust 
emission seen in an SNR. WISE data at 22 /im do not further constrain the models discussed later, 
and 12 p,m data are sensitive to a poorly-understood hot dust component, likely from very small 
dust grains which are stochastically heated. The IR model above predicts a flux of 100 mJy at 12 
fim, while the measured flux from WISE in this region is 430 mJy (albeit with a ~ 40% error). A 
similar excess of short-wavelength IR emission above a model fit to 24 and 70 /im was found in 



our Spitzer analysis of Kepler's SNR (IBlair et al. 



2007b . Emission at this wavelength would likely 



come from very small dust grains, and the physics of c ollisional heating and sputtering of such 



small grains is not well understood. See Section 3.3 of IBlair et al. 



(2007) for further discussion 



of the difficulties in modeling short wavelength IR emission. However, this emission from very 
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small grains should have little effect on the overall properties derived from the IR fits. If the 70 
fim data were ignored, and only the 12 and 24 ftm data were fit with a single-temperature grain 
model, the temperature of the grains would be 190 K, and the radiating mass would be only 2.4 

5 x 10" 7 M Q . Nevertheless, it is worth noting that significant work still needs to be done to 
understand short-wavelength IR emission from dust in SNRs. 



3.6. Dust-to-Gas Mass Ratio 



Dwek et al.l (|1987l) provide theoretical predictions for the ratio of the total IR luminosity to 
the 0.2-4 keV X-ray luminosity, a quantity that they refer to as the IRX ratio. At the temperature 
of the NW region, this ratio should be ~ 35. Our dust model calculates the amount of sputtering 
of grains in the post-shock gas, and for the NW region we find that 19% of the total dust mass 
has been destroyed. Accounting for this mass, and considering only the unabsorbed thermal 
component (vpshock) of our X-ray spectral fit described above, we find an IRX ratio of 1.65, 
lower than the predicted value by a factor of 20, which implies a dust-to-gas mass ra tio for the 



(119871) . The IRX 



pre- shock ISM lower than the standard Galactic value assumed by iDwek et al. 
ratio for the SW region is 1.4. 

We can approach this issue in another way. The X-ray EM for the thermal component in the 
NW region is 2.15 d| 5 cm" 3 M Q . Since we have (from the IR fits) an independent determination 
of the post-shock density, we can divide this out of the EM and get the amount of gas that has 
been swept up by the shock, which is roughly a solar mass. Normalizations to the IR spectra give 
the radiating dust mass, and dividing these two gives another measure (also distance-independent) 
of the dust-to-gas mass ratio, 1.37 x 10" 3 , where we have accounted for 19% of the mass in grains 
being destroy ed via sputtering. The typica l value for the dust-to-gas mass ratio in the Milky Way 
is 7.5 x 10" 3 (|Weingartner & Draindl2001|) . We thus confirm our finding from above of a low 
IRX ratio, with a dust-to-gas mass ratio that is roughly a factor of 5 too low. The discrepancy 
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between an IRX ratio that is 20 times too lo w and a dust -to-gas mass ratio that is 5 times too 



low primarily results from the assumption by 



Dwek et al. 



(I1987|) that the plasma is in collisional 



ionization equilibrium (CIE). In models for the X-ray emission from RCW 86, the plasma is in a 
non-equilibrium ionization state, which significantly raises the amount of X-ray emission when 
compared to a CIE model. For example, if kT = 0.35 keV and r, = 7 x 10 10 cm" 3 s, a plane-shock 
model will produce ~ 3 times more X-ray emission than a Raymond-Smith plasma in CIE. 

We arrive at a similar conclusion in the SW, where the dust/gas mass ratio is found to be 1.3 
x 10" 3 , but note that both of these methods will miss any dust that is too cold to emit at 70 jum 



Simi 



2007 



2006 



ar dust/gas mass ratio shortfa^ 



Arendt et al. 



2011). 



20101 : 



Lee et al. 



s have been observed in SN Rs in both the G alaxy (Blair et al 



2009) and the LMC (Bo rkowski et al. 



2006 



Williams et al. 



4. Progenitor Type of RCW 86 



The nature of the progenitor SN of RCW 86 is uncer tain, with most a uthors favoring a CC 



Westerlundl (119691) . who placed the 



origin. Support for the CC model primarily comes from 
location of the SNR within an association of ten B-type stars at a distance of 2.5 kpc (although 
one of these stars is likely a foreground object, and another has a doubtful association based on 
its distance). These stars are spread out over a much larger region of space than the periphery 
of the SNR. Recent optical studies of the remnant have used line kinematics of the narrow Ha 



components of various shocks in the r emnant, along with Ga 



distance to the remnant at 2.3-2.8 kpc (|Sollerman et al. 



2003; 



actic rotation curve s, to place the 



Rosado et al. 



1996), confirming that 



the remnant and the OB association are at approximately the same distance. 

However, this proximity to the association does not preclude a Type la progenitor much older 
and less massive than the B-type stars found in this region of our Galaxy. This is illustrated by 
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N103B, a young remnant in the Large Magellanic Clou d known to be t he res ult of a Type la SN, 



that is apparently associated with the H II region N103. 



Badenes et al. 



(2009) examined the local 



star formation history near N103B and found a prominent extended peak between 100 and 50 
Myr. A mixture of young and old stellar populations is likely present near RCW 86 as well, so 
both Type la and CC explosions may occur at this particular location within our Galaxy. 



Secondary support for the CC origin comes from the suggestion, first made by 



Vink et al. 



(|1997|) . that RCW 86 is an explosion into a cavity. Cavity explosions are typically associated 



with CC SNe and the strong winds that some m assive stars may blo w off prior to explosion 



However, as shown in 



Hachisu et al 



(1996) and 



Badenes et al. 



(|2007h . accreting WD systems, the 



progenitors of single-degenerate Type la SNe, can also create substantial wind-blown cavities. 

Both the arguments in favor of a CC origin are valid, but are circumstantial and based on 
general associations rather than specific observations of the remnant. Here, we critically examine 
available observations of RCW 86, and find that the remnant is more likely the remnant of a Type 
la SN. 

1. B aimer-dominated shocks. The entire periphery of the remn ant shows H a emission from 



non-radiative shocks encountering at least partially neutra l material (|Smith 



for this Ha emission is described in 



Chevalier et al. 



1997). The mechanism 



(1980). At the densities derived above from 



IR observations, H recombination times are ~ 10 5 yr. Most CC SNe ionize the surrounding 
medium in one of two ways, either from the radiation from the shock breakout (in the case of a 
red supergiant) or before the SN from the ionizing flux from the progenitor star itself (in the case 
of a Wolf-Rayet star). A hot, compact progenitor of a stripped CC Type Ib/c SN would produce 
relatively few ionizing photons in the explosion, but would have a strong ionizing flux before the 
explosion, as the progenitors of these SNe are massive O-type stars during their main-sequence 
lifetime, which proceed through a luminous blue variable (LBV) phase to shed their H (and 
possibly He) envelope, becoming Wolf-Rayet stars. Less massive (< 25 M & ), cooler stars will 
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to ionize material 



Chevalid(2005) 



be red supergiants (RSG) at the time of the explosion. While they are too coo 
beforehand, their large size results in a large ionizing flux at shock breakout, 
discusses the ionizing radiation from the shock breakout of a CC SN. For the amount of mass that 
the shock has encountered in RCW 86 (as constrained from model fits to IR data; see Section [5]), 
explosions ofRSGs ( R ~ 600i? o ), including Type IIP, Type IIL, and even Type lib SNe like 1993J 



(Woo slev et al 



19941) are likely ruled out, as they produce more than enough ionizing radiation. 



The progenitors of Type Iln SNe are general 



stars which evolve through an LBV phase (IGal-Yam et al 



ly thought to be extre mely mas sive (> 80 Mq) O-t ype 



20071) . although 



Smith et al 



(120091) 



point out that LBVs may not be the only possible progenitors of Type Iln SNe, and that they may 
arise from red supergiants as well. In either case, the mass contained in the CSM is enormous, 
and there is no evidence in RCW 86 for interaction of the blast wave with a dense CSM. The only 
remaining CC scenario which cannot be ruled out from this perspective is an explosion of a blue 
supergiant (BSG) progenitor with a massive H envelope, such as SN 1987A, producing only a 
modest number of ionizing photons, but see points 2 and 3 below. 

2. Fe abundance in ejecta. Even if some CC models cannot be r uled out from argument 1, all 



1995). 



CC SNe produce large amounts of O with respect to Fe in the ejecta (IWoosley & Weaveil . 
There is no evidence for super-solar O abundances anywhere in the remnant; both the spectra 
from the NW and the SW are consistent with typical ISM abundances for O. Fe L-shell emission 
is consistent with solar abundances as well (see TableQ]), but Fe Ka emission cannot be explained 
by this model. The properties of the fits listed in Table \T\ indicate that we are seeing shocked 
ISM. The 6.4 keV line observed in both regions is mos t easily explained b y a high-temperature 
(5 keV), low Tf (~ 1 x 10 9 cm" 3 s) plasma. As R02 and lUeno et al.l (|2007|) point out, a plasma at 
this temperature requires collisionless heating of electrons at the reverse shock. An O-rich and 
Fe-rich plasma in such a low ionization state does not produce L-shell emission from Fe, but H 
and He-like O K-lines are expected to be prominent, even if the ionization age of the plasma is 
lowered to a few x 10 8 cm" 3 s. A study of X-ray emission in the SW using Suzaku found a lower 
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limit of 0.05 d\ 5 M of Fe from a region cov ering 7% of the remnant's surface area (lUeno et al. 



2007|) . A similar study of the eastern limb by 



Yamaguchi et al. 



d2008al) found > 0.06 d\ 5 M & of 



Fe (with a possible upper limit an order of magnitude higher) from a region covering 2.5% of the 
remnant. As we show in Section[5l we find 3 x 10" 3 M Q in the NW, from a region containing 1.5% 
of the remnant. 

Can oxygen simply be hidden? This would require a very low ionization state of the gas. 



Ueno etal 



(2007) derive their Fe mass listed above by assuming that the electron density in the 
shocked ejecta is 0.1 cm" 3 , and that the ejecta were shocked 1000 years ago. They consider this 
electron density to be an upper limit, as anything higher would raise the ionization timescale of 
the gas to levels (> 3 x 10 9 cm" 3 s) that would result in strong Fe L- shell lines. Lowering the 
density to the point where oxygen lines would be weak, say, at r t = 10 8 would raise the amount 
of Fe in the SW region by a factor of 30, to 1.5 M . While this amount of Fe appears excessive 



even for a Type 



Yamaguchi et al 



a SN, it would also rule out all but the most exotic CC explosions. In the NE, 
(2008a) fit the centroid of the Fe Ka line with Suzaku, and found r, = 2.3 x 10 9 



cm 3 s, a sufficient timescale to produce strong O lines, yet none are reported. 



3. BSG Incompatibility with Observations. 



Gvaramadze & Vikhlininl (|2003|) proposed that 



the remnant is the result of the explosion of a BSG, moving through a bubble that the progenitor 
created during its main sequence lifetime. A large amount of mass may have been expelled at the 
RSG stage, facilitating transformation into the BSG. While the large amount of Fe seems to rule 
out any CC progenitors, there are two other arguments against this particular scenario. First, a 
slow wi nd blown at the RSG stage is expected to be nitrogen-rich, as obs erved in the CSM of SN 



1987A ( McCrav 



1993), but no nitrogen enhancement has been reported (|Leibowitz & Danziger 



19831) . Secondly, the explosion of a massive blue supergiant, su ch as the proge nitor of SN 198 7A, 



eads to the formation of an iron-nickel bubble (|Li et al. 



1993). As shown in 



Williams et al. 



(2008J), an iron-nickel bubble forms dense clumps of O-rich material, which should show strong 
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optical or X-ray emission. Again, no O-rich ejecta have ever been reported in this remnant. The 
failure of this elaborate CC scenario demonstrates the difficulty in finding a plausible CC model 
for RCW 86. 

4. Lack of central po int source or pulsar-wind nebula. While CC explosions, such as 



Cas A (IPavlov et al 
not. 



Fesenetal 



2000), often leave behind a compact remnant, Type la explosions should 



1979Q conducted an optical survey of several known historical SNRs in the 



Galaxy, finding no central point source in RCW 86 consistent with a compact remnant down to a 
magnitude of 18.5. iGvaramadze & Vikhlininl (|2003|) report the possibility of a point source from 



a compact object in the Chandra observations of the SW shell. However, this source is located 
near the edge of the remnant, and if associated with the supernova of 185 A.D., the transverse 
velocity of the sour ce would be in excess of 5000 km s" 1 , even assuming an off-center explosion. 

(|20041) conducted a search with Chandra of the interior of the remnant, finding no 



Kap lan et al. 



neutron star (down to a limit of one-tenth the luminosity of the Cas A neutron star) or pulsar-wind 
nebula. They assume the maximum reported distance of 2.8 kpc, and note that their limits would 
be even more constraining if the remnant were closer. However, their Chandra field of view does 
not cover the entire remnant, and they are unable to rule out less likely compact obj ects, such as 



massi ve neutron stars, neutron stars with exotic interiors, and quiescent black holes. 



Kaspi et al. 



(|1996h find no radio pulsar down to a limit of 1.3 mJy at 436 MHz and 0.2 mJy at 1520 MHz. No 



pulsar- wind nebula has ever been reported in the radio. 

Despite the previously mentioned proximity to an OB association and the suggestion that 
RCW 86 is the result of a cavity explosion, the specific characteristics of the remnant are more 
easily explained by a Type la explosion. We cannot definitively rule out all CC progenitors, but we 
can say that any such event would be a very unusual CC SN. It would have to produce significant 
amounts of Fe relative to O in the ejecta, while being a weak source of ionizing radiation, either 
before the SN or during the shock breakout phase. Its pre-supernova winds would be inconsistent 
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with the nitrogen-rich wind expected in the transition from an RSG to a BSG, and any compact 
remnant left behind would have to be less than one-tenth as luminous as the neutron star in Cas A. 

It is possible to confirm or refute the conclusion of RCW 86 being a Type la SN. Light 
echoes from historical SNe can reflect the spectrum of the SN event itself, and matching of 
detec ted echoes with SN spectral templates has been used in several cases to determine the SN 



type (IKrause et al 



2008). Detection of a light echo spectrum from this remnant that could be 



matched to a Type la template would confirm our conclusion. A caveat to this method is that 
not all thermonuclear explosions of white dwarfs look like Type la SN. A notable case of this 



is S N 2002ic, w 



SN (IWang et al. 



lerea strong H emission line was detected in the spectrum of a thermonuclear 



2004). Alternatively, detection of a compact remnant (neutron star, pulsar-wind 
nebula, or black hole) that is definitively associated with the remnant would clearly imply a CC 
explosion. Finally, although the remnant does not belong to the "O-rich" class of SNRs (a group 
of approximately half a dozen SNRs categorized by strong optical oxygen emission), it is still 
possible that deep X-ray observations of the center of the remnant could show oxygen in the 



ejecta. 



5. Hydrodynamic Modeling 



Perhaps the most puzzling aspect of RCW 86 is its large size, given its likely association 
with the supernova of 185 A.D. At a distance of 2.5 kpc, the diameter of the remnant is ~ 25 pc, 
requiring an average shock speed of > 7000 km s" 1 . Measured shock speeds vary greatly within 
the remnant. In the regions of interest considered in this paper (where the sh ell is detected at bot h 



24 and 70 /mi), s 



Ghavamian et al. 



lock speeds from optical spectroscopy are 500-700 km s 1 (|Rosado et al. 



1996; 



200 ll) , while shock speeds in the eastern limb from th e proper motion of X-ray 



emitting filaments have been reported as high as 6000 (±2800) km s 1 (IHelder et al. 



2009). The 



detection of thermal X-ray emission at ionization timescales of 6 - 10 x 10 cm s distinguishes 
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this remnant from SN 1006, another young remnant of a Type la explosion that is very large for 



its age . There, ionization timescales are shorter by nearly an order of magnitude dYamaguchi et al 



2008b) 



Badenes et al.l (|2007l) considered models of Type la explosions into bubbles blown by 
accretion wind outflows from the progenitor system, noting that most known Type la Galactic 
SNRs are inconsistent with this scenario. However, they note that RCW 86 may be the exception 



to this, and could in fact be the resu 



the winds modeled by Badenes et al. 



t of ju st such an event. The mechanism for the creation of 



Hachisu et al. 



(I1996h . A high mass 



(|2007h is described in 

accretion rate from a red-giant (RG) or main-sequence companion onto a white dwarf (WD) can 
lead to the formation of a large, loosely bound envelope surrounding the WD. This envelope 
is optically thick, and continuum-driven winds from the WD can eject several tenths of a solar 
mass from the system. The wind stops when the mass transfer rate falls below ~6x 10 -7 M Q 
yr" 1 , which can be as long as 10 5 - 10 6 yr before the WD accretes enough mass to reach the 
Chandrasekhar limit. 

We perform hydrodynamic simulations of Type la explosions into various surroundings 



using VH-1, a numerical hydrodyna mics code developed by J. Bl ondin and collaborators based 

"I 



on the piecewise parabolic method of Cole 



density profile of 



Dwarkadas & Chevalier ( 



la & Wo odward! (|1984|) . We use the exponential ejecta 



1998|) . and assume (except where otherwise noted) an 



explosion energy of 10 51 ergs and an ejected mass of 1.4 M Q . Densities reported below are in units 
of either H atoms cm" 3 , where we use a scaling of n Q = po/2.34 x 10" 24 gm cm" 3 , appropriate for 
an ISM with a 10: 1 H:He ratio, or in amu cm" 3 . We seek solutions which approximately reproduce 
the known age, radius, and velocity. Radius is a difficult quantity to determine for the somewhat 
unusual morphology of RCW 86. We chose a point at the approximate geometric center of the 
remnant and measure (assuming D = 2.5 kpc) radii of 12 and 13 pc to the NW and SW filaments, 
respectively. Since the NW represents the simpler geometry of the two filaments, we chose 12 pc 
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as the radius to be reproduced in the hydro simulations. 



5.1. Explosion into a Uniform Medium 

We beg in with the simplest case, an e xplosion into a medium with constant density, as 
modeled in Dwarkadas & Chevalien (|1998f) . In a 1-D hydro model, assuming the standard 
explosion parameters above, the forward shock can reach the desired radius in 1825 years by 
encountering a constant density of n = 0.033 cm" 3 . In addition to the density in this model being 
more than an order of magnitude below that required to produce IR emission in the NW and SW, 
the current forward shock velocity in this model is v s = 3000 km s" 1 , well above that measured 
from optical spectroscopy. 

An alternative model can reach 12 pc with a current shock velocity of 600 km s" 1 by 
encountering a medium of n = 0.75 cm" 3 . While this model is consistent with both the measured 
shock velocities and the inferred densities from IR, the time required for this shock to reach 
the observed radius is 7000 years, ruling out an association with SN 185 A.D. A third option 
is to require the shock to decelerate to 600 km s" 1 in 1800 years, but this model requires an 
ambient density of n = 40 cm" 3 and places the current radius of the forward shock at 3.2 pc. The 
distance to the remnant derived from this value would be 650 pc, well inside any previous distance 
estimates. Similar calculations show that reasonable under- or over-energetic SNe likewise fail to 
reproduce anything close to the observed remnant. 



5.2. Wind-Blown Bubble Model 



Given the results of models discussed above, we can rule out the case of a blast wave from an 
explosion in 185 A.D. encountering a uniform ambient density. The next level of complexity that 
we consider is a low-density, uniform bubble (region I) surrounded by higher-density, uniform 
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ISM (region II). A constant wind blown from the progenitor system would lead to an r" 2 density 
profile, but a uniform bubble could be created by simply switc hing off the wind a fe w thousand 
years prior to the explosion; delay times prior to explosion in iBadenes et al.l (|2007|) are well 
in excess of this in most cases. Since our IR modeling gives us a direct measurement of the 
post-shock density of the material recently shocked, we have one more constraint to use in these 
models, and we fix the ambient density in the ISM to the post-shock density from IR fits (~ 2 
cm" 3 ) divided by a compression ratio of 4, i.e. no,// = 0.5 cm" 3 . With this number fixed, the 
number of free parameters in our bubble model is only two: the density within the bubble, which 
we define as n j, and the location of the bubble "wall," which we define as R. 

We find a good fit to the observed radius and age from a model with n ,/ = 0.002 cm" 3 and R 
= 10.8 pc. The current shock speed in this model is 740 km s" 1 , which roughly agrees with the 
600-700 km s" 1 shock speeds measured from optical data in the NW and SW, respectively. We 
show the density profile from a one-dimensional (1-D) simulation of this scenario in Figure |6] 
In this model, the shock races through the low-density bubble in 725 years, reaching the shell 
with a velocity of 8900 km s" 1 . Upon hitting the shell o f density contrast 2 50, the transmitted 



Borkowski et al. 



(|1997|) provide an analytic 



shock speed into the dense material is 1300 km s" 1 . 
estimate for the transmitted shock speed when transitioning from a bubble to a dense shell. If 
5 is the density ratio, no,///"o,/, between the dense shell and the bubble, then the pressure jump, 

P = Pi/Pi, is 



P = 6(1 + 1. 701<T 1/2 -0.40185" 1 + 0.2274<T 3/2 -0.08745" 2 )" 



(1) 



and the transmitted shock speed is 



V t = ((3/5)^ 2 V B 



(2) 
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where V B is the speed of the shock when it reaches the wall. With 5 = 250 and Vb = 8900 km s 1 , 
the transmitted shock speed in this approximation is 1250 km s" 1 . 

The shock has been in the dense ISM for 1 175 years, sweeping up a total mass of ~ 35 M & . 
The reverse shock, meanwhile, has shocked virtually all of the ejecta (1.2 M ). The ionization 
timescale of the forward shocked material implied by this model, r, = 8.9 x 10 10 cm" 3 s, agrees 
quite well with the values of 9.8 x 10 10 cm" 3 s and 6.0 x 10 10 cm" 3 s derived from X-ray fits to the 
NW and SW. We can use a semi-analytic approximation to check the final shock speed (740 km 
s" 1 ), by assuming that the internal energy in the shocked ejecta is the dominant energy component, 
significantly greater than the internal energy in the shocked shell or the kinetic energies of the 
ejecta and shell material. In this approximation, the forward shock is completely driven by the 
pressure of the shocked ejecta, and the shock speed is given by V s = (Ap/3p s ) l l 2 , where p is the 
thermal pressure, p = E/2ttR 3 . For an explosion energy of 10 51 ergs and R = 10.8 pc, p = 4.3 
x 10" 9 ergs cm" 3 . In this approximation, Vs = 700 km s" 1 . In reality, the internal energy of the 
shocked ejecta, obtained by numerically integrating the output of the hydro model, is ~ 40% 
greater than the sum of the kinetic and internal energies of the shocked ISM, so this approximation 
is not entirely valid, but the general agreement is encouraging. 

Once the assumption of a uniform ambient density is abandoned, the parameter space for 
possible density structures that the shock can encounter is vast. The solution above of a shock 
running through a low-density bubble and encountering the much higher density ISM is not 
unique, but it is not our aim to explore every possible scenario. This model represents the simplest 
solution consistent with the available observed c onstraints. It is interesting to note that while 



we do not use the models of 



Badenes et al. 



(2007J) as a constraint in any way, the best-fit model 



listed above looks remarkably like model HI from their paper, which consists of a white dwarf 
of 1 M Q in a binary system with a donor star of 2 M , with an initial period for the system of 
2 days. The total mass lost by the system is 0.15 M , with an outflow velocity of 1000 km s" 1 , 
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and the time between the onset of the final mass transfer and the SN explosion is 5 x 10 5 yr. The 
"bubble" created by this system has a roughly constant density of 0.001 cm" 3 with a shell wall 
located at ~ 17 pc. The density in the ISM is fixed at 0.42 cm" 3 , with the exception of a shell at 
the bubble-ISM interface with a density higher by a factor of ~ 2. We assume that the blast wave 
is contained within this shell. 

This model also provides an explanation for the Fe Ka line seen in the X-ray spectrum from 
the narrow extraction region in NW, where the elliptical region is quite narrow compared to the 
shock radius. The extraction region, shown in Figured! begins at a radial distance of 93% of the 
forward shock radius. Shocked ejecta cannot exist this far out in the 1-D model, but this is only 
indicative of the limitations of modeling a multi-dimensional system in 1-D. In two dimensions 
(2-D), the instabilities generated by the interaction of the blast wave with the density profile in the 
bubble model allow the ejecta to pile up quite close to the forward shock. We show the output of 
a 2-D model at the current age of the remnant in Figure H which shows significant amounts of 
shocked ejecta at a radius of > 93% of the forward shock radius. 

We can also use the output of the hydrodynamic modeling, combined with X-ray 
spectroscopy, to estimate the amount of reverse- shocked Fe present in the NW extraction region. 
The EM, as fit to the Fe Ka line in the XMM-Newton spectrum for the Fe-only component, is 7.2 
x 10" 6 cm" 3 M Q . The hydro models give the density in the reverse shocked ejecta, which over the 
small range of interest (> 93% of the forward shock radius) is roughly constant at 1.38 x 10" 2 
amu cm" 3 . If this is all Fe, then n Fe = 2.5 x 10" 4 cm" 3 . If we assume 10 free electrons per Fe atom 
(to avoid ionizing Fe to the point of producing Fe L-shell emission), then an electron density of 
2.5 x 10" 3 cm" 3 leads to an Fe mass of 3 x 10" 3 d\ 5 M Q , a lower limit given the approximations 
used. This is a relatively small amount of iron, but the size of the extraction region for the NW 
covers only abo ut 1.5% of the s urface a rea of the remnant. As we mention in Section 51 similar 



calculations by 



Ueno et al 



(2007J) and 



Yamaguchi et al 



d2008ah find a lower limit of 0. 1 1 dj 5 M Q 



of Fe from less than 10% of the area of the remnant. 



5.3. Off-Center Explosion 

Despite the general agreements listed above between our 2-D hydrodynamic model and the 
measured properties of the remnant, one thing this model cannot account for is the extremely high 
shoc k speed reported f or the NE limb, found from Chandra proper motions to be 6000 ±2800 km 



s 1 (Helder et al. 



3) 



2009). Our hydro model assumes spherical symmetry for the distance from the 



explosion to the shell wall, and this symmetry must be broken to explain shock speeds that differ 
by an order of magnitude from one side of the remnant to the other. It is beyond the scope of this 
paper to provide a detailed multi-dimensional model to explain the morphological "squareness" 
and varying shock speeds of RCW 86; such a model would undoubtedly be highly complex. 
Rather, we seek the simplest case that can broadly explain the observations, that of an off-center 
explosion into a wind-blown bubble, in which the shock has not yet reached the bubble wall in the 
NE. 

We assume that the wind-blown bubble is symmetric with a radius of 12 pc. We then vary the 
distance, R , from the center of the bubble to the center of explosion. R is the only free parameter 
in this model; we assume the same bubble and shell densities as above (n Q I = 0.002 cm" 3 and n 0H 
= 0.5 cm" 3 ). In this model, the distance, D, from the center of explosion to any point on the wall 
of the shell is given by 



D(R J) = (Rl-2R Rcos 6+R 2 )^ (3) 

where R is the fixed radius of the shell wall (in this case, 12 pc) and 6 is the polar angle between 
the center of explosion and a point on the shell wall. A diagram of the geometry of this model is 
shown in Figured In 2-D, we model this as a semi-circle from to ir radians, using a spherical 
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geometry with reflecting boundary conditions at and ir. 

Within this symmetric shell model, a significant offset of the explosion from the center of the 
bubble is required for a shock that roughly reproduces the observables (shock velocity, density, 
age) in the SW, yet still allows for the extremely high shock speeds in the NE. To avoid hitting 
the far side (9 = ir) of the shell wall in 1825 yr, R must be at least 7 pc. Assuming this value, the 
shock encounters the shell wall on the near side (6 = 0; the SW corner of the remnant) in ~ 300 
yr, and reaches a final radius (with respect to the center of the remnant) of 13.5 p c with a veloc ity 



Rosado et al 



of 1 000 km s 1 . The hydrodynamic output of this model is shown in Figure [9l 
(|1996|) estimated the shock velocity in the SW region to be 800 km s" 1 . At 9 = ir/2, which would 



roughly correspond to the location of the NW filament, the shock hits the filament in ~ 1000 yr, 
reaching a final radius (with respect to the center of the remnant) of 12.7 pc with a velocity of 550 
km s" 1 , numbers which agree quite well with observations. This model implies that the shock has 
only recently encountered th e shell in the E and S E. The ionization timescale of the shocked ISM 



in the E, as reported by 



Yamaguchi et al 



(2008a), was 7.7 x 10 u cm J s. This is rather high for 



a recently shocked shell, but the authors point out that densities in this region are likely much 
higher than the rest of the shell, and that the ionization timescale of the reverse shocked Fe implies 
a shock age of < 380 yr. Thus, there is a broad consistency of observable parameters with the 
off-center explosion model, including the transition from thermal to nonthermal X-rays in the NE 
corner of the remnant, where, in our model, the shock has not yet hit the bubble wall. Getting the 
progenitor system away from the center of the bubble does not require excessively high velocities. 
If the delay time between the onset of the fast wind and the SN explosion is as little as 10 5 years, 
a velocity of 65 km s" 1 would cause an offset of 7 pc. Delay times of 10 6 years or more require 
velocities of only a few km s" 1 . 
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6. Synchrotron X-rays 



The nonthermal X-ray emission from RCW 86, especially in the SW, differs in some 
important ways from that seen in other young SNRs. Most young SNRs show synchrotron X-rays 
concentrated in "thin rims" at the periphery (see Reynolds 2008 for a review). These rims are 
commonly assumed to be limited in thickness by the distance an electron can advect with the 
post-shock flow in a synchrotron loss time (Vink & Laming 2003; Parizot et al. 2006). However, 
in RCW 86, most synchrotron emission is more diffuse, except for the E and NE edges. The 
spectrum in the SW can be well described by emission from a power-law electron distribution 
with exponential cutoff (model srcut in XSPEC) (R02). The characteristic "rolloff" frequency v roa 
(critical, not peak, frequency of electrons with the e-folding cutoff energy E m ) is about 10 17 Hz 
(R02), implying E m = A0{B w )~ l l 2 TeV, where B i0 =5/10 fiG. R02 showed that, given a geometric 
effect that can speed up particle acceleration in perpendicular shocks (Jokipii 1987), shocks of 
the speed of the optical nonradiative shocks in the SW (~ 600 km s" 1 ) could accelerate electrons 
to this energy. However, without this obliquity effect, shock velocities in excess of 1000 km s" 1 
are more typically required. Such shocks have not been present in the SW since the collision 
about 1500 yr ago of the blast wave with the bubble wall. However, at that time the blast-wave 
speed was of order 10,000 km s" 1 in our scenario, and acceleration to very high energies (100 TeV 
or above) would have been possible due to the strong dependence of maximum energy on shock 
speed (^max oc V v 2 for loss-limited acceleration). In this case, presumably a "thin rim" morphology 
would have resulted, but the rolloff frequency would have been considerably higher than observed 
today. But electrons with energies below E m of that time have longer lifetimes: the time for an 
electron of energy E TeV to lose half its energy is ty 2 = 1.3 x lO 5 .^ -1 !^ Y r - Combining this with 
the relation between energy and characteristic frequency, we find that electrons with characteristic 
frequency v mVi can survive for a time t if they radiate in a magnetic field less than 
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(We note that if B < 3.3 /iG, the magnetic-field strength with energy density of the cosmic 
microwave background (CMB), then inverse-Compton losses from CMB photons will dominate 
synchrotron losses. However, this is unlikely to be the case in RCW 86 or other young SNRs.) For 
RCW 86, electrons with v mU ~ 10 17 Hz can survive for 1500 yr if B < 22 /iG. 

Synchrotron emission from electrons accelerated over 1500 years ago would not be expected 
to be concentrated in thin rims; it could occupy the entire region of shocked bubble material, 
and would appear in projection as diffuse emission extending back into the remnant interior, as 
observed. For younger remnants like Tycho, thin rims (according to the synchrotron-loss scenario) 
imply much higher magnetic fields which would deplete all electrons capable of producing 
X-rays. However, in RCW 86 and perhaps other cavity explosions, it is natural to expect diffuse, 
long-lasting synchrotron X-ray emission, as long as magnetic field strengths remain relatively 
low. Acceleration to very high energies, while the blast wave is still inside (or within) the bubble, 
would produce a population of electrons that could go on emitting a synchrotron spectrum up 
to X-ray energies, even if the subsequent shock transmitted into the dense surroundings rapidly 
becomes too slow to accelerate electrons to TeV energies. 

A detailed analysis of the X-ray synchrotron emission from all around the periphery of RCW 
86 is beyond the scope of this paper, but should be undertaken to test this suggestion and to clarify 
the current state and past history of electron acceleration in RCW 86. 

7. Gamma-Ray Emission 

We can use our IR fitted densities to make an order of magnitude estimate regarding the 
gamma-ray emission observed. In the hadronic model of gamma-ray emission, cosmic-rays 
escaping the shock collide with thermal gas particles, producing tt° particles which decay into 
gamma-rays. The 1-10 TeV spectrum from H.E.S.S. is fit with a photon index of T = 2.41 
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(lAharonian et al 



2009(). If extrapolated back to 100 MeV, this gives a 



l umino sity, L 7 j om<>v = 1.19 



Drurv et al. 



(I1994|) . the energy density 



x 10 39 photons s" 1 . Using the gamma-ray emissivity from 
in relativistic particles for a spherical SNR required to produce a 7-ray luminosity of L 7 by this 
process is 



_ ~ n _i7 / L l W0MeV \ f d \ 3 / 9 rad \ 3 / n H \ 1 _j _ 3 

where L 7 jooMev is the luminosity, in photons s" 1 , of 100 MeV 7-rays, d is the distance to the 
remnant, 9 rad is the radius of the remnant in arcminutes, n H is the post-shock proton density, and 
q 7 is the production rate of 7-rays, in units of s" 1 erg" 1 cm 3 H" 1 . We take a representative value 



from 



Drury et al. 



19941) of q 7 = 5 x 10 -14 and calculate u rel = 5.1 x 10 -8 ergs cm" 3 for RCW 86. If 
we make the (admittedly crude) assumption that the NW filament typifies the entire remnant, then 
the total energy density (pov 2 ) of the thermal gas is 4.2 x 10" 9 ergs cm" 3 . In this case, u re i/p v 2 s = 
12, requiring an order of magnitude more energy going into relativistic particles than is available 
in the shock. A more realistic value for this ratio, 0.1, requires a flatte ning of the gamma-ray 



Aharonian et al. 



(2009), 



spectrum below 1 TeV to T = 1.9. We thus confirm the conclusions of . 
who based their argument for the low-energy hardening of the 7-ray spectrum on the total energy 
from the SN of 10 51 ergs. 

Any particle acceleration occuring in the remnant would rob the post-shock gas of energy, 
raising the shock compression ratio across the shock. We cannot determine this ratio directly, 
since we do not have a direct measure of the emitting volume of the NW region. However, since 
we have the post-shock density and total shocked gas mass, the volume for a "standard" shock 
with a compression ratio of 4 would be 1.52 x 10 57 cm 3 . Assuming D = 2.5 kpc, the projected 
radii of the elliptical region are 4.7 and 0.62 pc for the major and minor axes, and the line-of-sight 
depth through the filament would be 4.3 pc, a factor of 7 greater than the radial width. This is the 
minimum line-of-sight, as anything shorter would imply a CR < 4. 
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8. Concluding Remarks 

Spitzer and WISE mid-IR images show the complete shell of RCW 86, with the dominant 
source of emission being warm dust grains heated and sputtered in the post-shock environment. 
A combined IR/X-ray analysis yields post-shock densities in the non-radiative regions of the NW 
and SW of n H = 2.0 cm" 3 and 2.4 cm" 3 , respectively, with ~ 20% of the dust in these regions 
destroyed via sputtering. Based on the IR/X-ray flux ratio, we find that the dust/gas mass ratio in 
the pre-shock medium is lower by a factor of ~ 5 from that typical for the Galaxy. The post-shock 
gas densities derived from IR observations place strong constraints on the total amount of mass 
that has been encountered by the shock, and favor a young age for the SNR, strengthening the 
case for an association with SN 185. 

X-ray and optical evidence points to a Type la origin for the progenitor of RCW 86. The 
presence of Balmer-dominated shocks around the periphery of a young SNR means that the 
blast wave is encountering neutral material, and most CC SNe ionize significant portions of the 
surrounding ISM, much more than is observed here. There are SNRs which are believed to be the 
result of explosions within a cavity and show Balmer emission at the shock front (e.g., the Cygnus 
Loop, believed to be the remnant of a CC SN), but these are older objects where the forward 
shocks have encountered significantly more mass (the swept mass of t he Cygnus Loop is like ly 



> IOOMq, based on the densities obtained by 



Sankrit et al. 



(12010 ) and 



Raymond et al. 



((20031)). 



Not all CC explosions can be ruled out because of the Ha emission, but the amount of Fe in the 
ejecta relative to O and other elements also indicates a Type la origin. The X-ray analyses of the 
amount of Fe relative to O in the ejecta have been taken from references already in the literature, 
but the constraints on the ionizing flux from the progenitor come from the new IR observations 
presented here. There is additional evidence to support the claim that RCW 86 is the result of 
a Type la SN, in that blue supergiant explosions, which may avoid ionizing the surrounding 
medium, are incompatible with optical observations, and that no compact remnant or pulsar-wind 
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nebula has ever been detected. 

Most intriguingly then, we find that if RCW 86 is the remnant of 185 A.D., and is also the 
result of a Type la SN, then the only self-consistent model to explain the IR, X-ray, and optical 
observations is that of an explosion into a cavity created by the progenitor system, a model 
which requires a single-degenerate progenitor. We used one and two-dimensional hydrodynamic 
models to arrive at this conclusion, finding that a simple model of an explosion in a low-density 
bubble surrounded by a higher density shell can approximately reproduce observed shock radii, 
velocities, and post-shock gas densities. The density contrast between the bubble and the shell is 
250. This model does not require an appeal to any special circumstances, such as an exceptionally 
recent encounter with the shell wall or a sub/super-energetic explosion. The high proper motion of 
the X-ray emitting filaments in the NE, and thus the inferred high shock speed there, requires the 
explosion to be off-center, so that the shock in that region has not yet encountered the dense shell. 
In this model, the blast wave began encountering the dense ISM ~ 1500 yr ago and continues to 
impact parts of the shell. 

The fast, synchrotron X-ray emitting shocks in the NE are still encountering ionized material 
within the bubble, and thus no Balmer emission should be expected from these shocks. In the 
places where the shock is encountering the shell of the bubble, it does so at an oblique angle. 
X-ray proper motion measurements are sensitive to the speed at which the shock moves as it 
encounters the shell, which is ~ 5500 km s" 1 in the 2-D hydrodynamic model. The broad Ha line 
width, on the other hand, will reflect the post-shock temperature within the dense shell, where the 
transmitted shock is considerably slower. In our hydro model, the shock speed of the Ha emitting 
shocks within the dense shell is ~ 800 km s" 1 . We did not adjust this model to match the values 
for the NE, nonet heless, these values are quite close to the 6000 and 1 100 km s" 1 shocks reported 



for that region by 



Helderetal 



(2009). 



Thus, this model offers an alternative explanation for the apparent discrepancy between the 
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shock speed meas ured by proper motions and that derived from optical line kinematics, which 



Helder et al. 



(|2009|) attributed to significant cosmic-ray accleration, with > 50% of the shock 
energy going into cosmic rays. While we do not rule it out, our model does not require any 
acceleration of cosmic rays at the shock front. 

The synchrotron emission in RCW 86 is complicated, with the SW and NW showing 
more diffuse emission than the NE, where rims are sharper. The bubble model offers a natural 
explanation for this. In the regions of the remnant where the shock encountered the shell long 
ago, the electrons were accelerated by much faster shocks than are currently present there. For 
reasonable values of the magnetic field, these electrons can survive to the current age of RCW 86. 

This idea may find applicability in another older SNR traditionally ascribed to a cavity 
explosion: the Cygnus Loop, where GeV emission has recently been detected by the Fermi 
satellite (Katagiri et al., submitted). Electrons accelerated early in the life of the Cygnus 
Loop, while its blast wave moved rapidly through the low-density bubble, could persist for 
the estimated age of order 10,000 yr (e.g., Levenson et al. 1998) with energies as large as 
10 B^l TeV. Such electrons could produce gamma rays by inverse-Compton scattering of various 
photon fields. As a rough example, 1 GeV photons could be produced from upscattering cosmic 
microwave background (CMB) photons (£~3x 10" 4 eV) by electrons with Lorentz factors 
7 ~ {hu out /hi/i n )~ 1 / 2 ~2x 10 6 , or electron energies of about 1 TeV. Electrons with these energies 
could easily survive from an early period of rapid particle acceleration as long as B < 25 fiG or 
so in the regions they inhabit. In a current magnetic field B m ~ 1, they would have characteristic 
frequencies of synchrotron radiation u c ~ 2 x 10 14 Hz, i.e., in the near-IR, so no nonthermal 
X-rays would be expected. 

Thus, RCW 86 may indeed be the first example of a Type la explosion in a system in 
which the progenitor carved a wind-blown bubble, and can provide new insight into the outflows 
from accreting systems. While the hydro simulations themselves cannot distinguish between 
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a Type la and CC progenitor system, they can rule out the case of a uniform ambient density. 



We find a general agreement with the parameter s of the bubb 



dynamics of the remnant and the simulations of 



Badenes et al 



e model necessary to explain the 



(|2007|) . who modeled the outflows 



of single-degenerate Type la progenitor systems. In such a system, the companion star would 
retain its generally high velocity after the explosion. A star in a binary system with a period of 
a few days can h ave an orbital velocity , V, of up to several hundred km s" 1 . As shown in the 



recent work by 



Kerzendorf et al. 



(|2009|) . high-resolution spectroscopy to determine the rotational 
period can test the viability of a potential companion star. The maximum angular distance that a 
companion star could travel from the center of explosion in 1825 yr is 47 V 3 oo d 2 .5 arcsec, where 
V300 is the velocity of the companion in units of 300 km s" 1 . Although the center of explosion is 
currently unknown, it should be possible to identify the companion star from a single-degenerate 
Type la SN. Another promising possibility would be to find the light echo of RCW 86 and obtain 
a spectrum of the SN itself. 
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Table 1 . X-ray Spectral Models 
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Note. — Errors quoted correspond to 90% 
confidence intervals. Abundances are from 



Wilms et al. 



(120001) . Models above do not in- 
clude Fe-only component to account for Fe K- 
a line at 6.4 keV. srcut parameters (a, v ro u- ff, 
& Fj ghz) are frozen to a fit performed only 
from 2-6 keV. Errors listed are for that fit 
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Table 2. IR Models 



IR Measurements and Modeling NW SW 
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Note. — Md MS ( is pre-shock dust mass encountered by blast wave, 
i.e., has already been corrected for destruction of grains via sputtering. 
Lir is total IR luminosity, from 1-400 fim. Shock age is calculated by 
dividing 77 by n e . 
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Fig. 1 . — Top Left: MIPS 24 /im image, with NW extraction region indicated as green ellipse; 
Top Right: MIPS 70 /im image; Bottom Right: XMM-Newton EPIC image, with 0.5-1.0 keV 
in red, 1 keV to 1.95 ke V in green, and 1.95-6.8 keV in blue; Bottom Left: NOAO/CTIO star- 



subtracted optical image (ISmith 



1997|) . with Ha in red and [S II] in green. Areas appearing yellow 



contain both non-radiative and radiative shocks. Cyan rectangle (lower left panel) marks Spitzer 
IRS spectral extraction region, as described in text. In all images, N is up and E is to the left. The 
remnant is approximately 40' in diameter. 





Fig. 2. — Top Left: IRAC 3-color image of SW region, with 8 /im in red, 5.8 /mi in green, 4.5 /mi 
in blue. Yellow regions highlight structures also seen at 24 /mi; Top Right: MIPS 24 /mi image 
with IRAC regions overlaid; Bottom Right: optical Ha and [S II] image, as in Figured! with region 
highlighting a purely non-radiative filament, used for IR analysis of the SW region, as discussed 
in the text; Bottom Left: EPIC-MOS2 X-ray image, with colors as in Figure [TJ 




Fig. 3. — Spitzer and WISE mosaic of the entire remnant, with Spitzer 24 /mi emission in red, 
WISE 12 /mi emission in green, and WISE 4.6 /mi emission in blue. The 22 /mi image from WISE, 
not shown, looks virtually identical to the Spitzer 24 /im image. We display the MIPS 24 /mi image 
here because of the sharper resolution and better signal-to-noise ratio. Only the brightest sections 
of the NW and SW are visible at 12 /im (note the slightly yellowish colors of these regions). 
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Energy (keV) 



Fig. 4. — EPIC-MOS 1 and 2 spectra, in black and red, respectively, of NW region, overlaid with 
a model of thermal and non-thermal emission, as described in Section [331 An additional thermal 
model to account for the Fe Ka line at 6.4 keV is included as well. Data are 3-cr binned for plotting 
purposes only. 
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Fig. 5.— EPIC-MOS image of the NW filament, with 0.3-0.7 keV emission in red, 0.7-1.3 keV 
emission in green, and 1.3-6.8 keV nonthermal emission shown in blue. The low-energy bands 
highlight thermal emission, while the high-energy band shows the location of non-thermal syn- 
chrotron radiation. Images have been smoothed with a 2-pixel Gaussian to highlight extended 
emission. 
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Fig. 6. — Density profile of the wind-blown bubble model, as described in the text, in one dimen- 
sion. Identified in the image are the forward shock (FS), the contact discontinuity separating the 
shocked ISM material from the shocked bubble material (CDl), the contact discontinuity separat- 
ing the shocked bubble material from the shocked ejecta (CD2), and the reverse shock (RS). The 
collision of the blast wave with a shell of density contrast 250 leads to density structures from 
transmitted and reflected shocks, which can also be seen in the figure. Pressure, normalized to 1 at 
CD2, is shown as a dashed line. 



Fig. 7. — Two-dimensional density structure in the wind-blown bubble model for an explosion 
located at the center of the bubble. Radius is normalized to 1 at the forward shock location (~ 12 
pc). Marked radii are forward shock (1) and reverse shock (0.63). Orange and red colors indicate 
shocked shell material, white represents shocked "bubble" material, and blues represent reverse- 
shocked ejecta. The Rayleigh-Taylor unstable contact discontinuity is traced by the border between 
white and blue. The collision with the dense shell brings the contact discontinuity much closer to 
the forward shock than is possible from a shock expanding into a uniform medium. Polar angle 
from 0.47T to 0.67T. 
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Fig. 8. — Geometry of off-center explosion model, as described in the text. The center of the 
spherically symmetric bubble is marked by a dot, and point of explosion is marked by a star. In 
this model, the radius from the center of the bubble to the inner wall of the shell is the approximate 
observed radius of the remnant, 12 pc. R is the distance from the center of the bubble to the center 
of the explosion, which is variable in this model. 
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Fig. 9. — Density structure in the off-center explosion model, as seen in the plane of the sky, shown 
at the current age of the remnant. The center of the explosion is marked with a black star, and the 
center of the bubble is marked with a yellow star. The distance between the two (R in the text) is 7 
pc. At 6 = 0, measured from the center of explosion, the foward shock has reached a radius of 13.5 
pc with respect to the center of the bubble, while at 6 = |, the shock is at 12.7 pc, again with respect 
to the center of the bubble. On the far side (9 = it), the forward shock is still in the low-density 
bubble, and has not yet reached the wall. Lack of an instability at the contact discontinuity is the 
result of the low angular resolution used in this model. 



